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The tunneling through an opaque barrier with a strong oscillating component is investigated. It is shown, that 
in the strong perturbations regime (in contrast to the weak one), higher perturbations rate does not necessarily 
improve the activation. In fact, in this regime two rival factors play a role, and as a consequence, this tunneling 
System behaves like a sensitive frequency-shifter device: for most incident partides' energies activation occurs and 
the partides are energetically elevated , while for specific energies activation is depressed and the transmission is 
very low. This effect is unique to the strong perturbation regime, and it is totally absent in the weak perturbation 
case. Moreover, it cannot be deduced even in the adiabatic regime. It is conjectured that this mechanism can 
be used as a frequency-dependent transistor, in which the device's transmission is governed by the external field 
frequency. 



When a quantum particle propagates through an opaque 
barrier in the classically forbidden domain, it tunnels. 
This conduct is manifested by the exponentially small 
transmission probability. Nevertheless, when some part 
of the barrier weakly oscillates the particle will be ac- 
tivated to higher energies, and the transmission will in- 
crease substantially ||ï|- However, when the temporal 
change is not merely a perturbation but rather a strong 
variation the dynamics are too complicated to be de- 
scribed by such a simple rule of thumb, and most of the 
interesting phenomena belong to this category. 

While it is well known that even a very weak external 
oscillating field may considerably increase the tunneling 
current |^,^ , changes which are comparable to the initial 
System parameters can cause elevator resonant activation 
and activation assisted tunneling charge quanti- 

zation pumping |Q , as well as wave function coUapse . 
That is, large changes reveal a wealth of physical phe- 
nomena. 

When the tunneling particle energy (Í7) is close to the 
potential barrier height (V), the dynamics become morè 
complicated since the perturbation amplitude (AT^) can 
exceed V — fi , i.e., the perturbation "strength" may be 
larger than the effective tunneling barrier. This case is 
extremely sensitive, since the dynamics are governed by 
two rival factors. On the one hand, the oscillation's am- 
plitude is so large, that for a finite segment of the oscilla- 
tion period the alternating potential blocks the partides' 
transmission. Hence, foUowing this reasoning, the tun- 
neling rate should decrease. On the other hand, energy 
quanta generated by the oscillating barrier can be ab- 
sorbed by the tunneling particle to assist it in the due 
course of tunneling. 



erel.g@kailight .com 



In this work, we investigate the interplay between these 
two rival factors, and show that the competition between 
them is responsible for the high sensitivity to the system's 
parameters. For example, it so happens that the incident 
partides are not always activated to higher energy: when 
their initial energy is equal to one of a series of specific 
energies the activation is frustrated, elevation to higher 
energies does not occur and the transmission is decreased 
accordingly. 

The problem of tunneling in the presence of an external 
time-dependent potential has been extensively studied 
||]. However, in this paper the strong perturbation 
regime is addressed and no approximate assumptions are 
made in the solution analysis. We give this case both 
an exact numerical solution, and an analytical solution, 
using an approximation, to show that the general conclu- 
sions can be deduced even in the slowly varying regime. 
However, we also show that the oscillations rate cannot 
be arbitrarily small, and in fact it should be larger than 
the spectral bandwidth of the resonances, otherwise (and, 
in particular, the adiabatic regime) the effect vanishes. 

The apparent underlying physics can be described as 
foUows: Due to the oscillations, the alternating barrier 
functions as an impenetrable barrier (or at least a very 
opaque one) only for a fraction of the oscillating period 
and therefore, a time-window (denoted tw) is formed in 
which the incident particle can propagate. Consequently, 
in order to prevent a destructive pattern, the incoming 
partides' frequency ,íí (in units where ?ï = 1) should be 
equal to an integer of both the oscillations' frequency lo 
and 2'k/t\y (i.e., Q — 27rmi/Tw and Q — for inte- 
ger mi and m2). Therefore, a selection rule for activation 
should be expected . However, this selection rule is rather 
complicated since tw is a function of the incident par- 
tides' frequency (íl). It turns out that the energies, for 
which activation occurs, have a nontrivial dependence on 
the integer nii (roughly like n\^^). 
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A few words should be added here about the plausi- 
bility of the efFect's practical implementation. The in- 
dustry has a special interest in resonant tunnehng de- 
vices(RTD), which allow miniaturizing electrònic circuits 
and improving their performance. Ordinary RTD are 
very sensitive to manufacturing processes, temperature, 
and impurities, and as a consequence a rehable device is 
highly costly. The presented tunnehng device aUows for 
the fabrication of a low-cost device where any resonant 
reíinement can be done by variations in the external field 
Jrequency, and no special geometry or manufacturing re- 
strictions are needed. 

In this paper we discuss the tunnehng dynamics of a 
beam of partides which are activated by strong harmònic 
perturbations. The tunnehng takes place through a very 
opaque potential barrier (high and wide). In order to 
have the barrier opaque at all times, we foUow Q and 
discuss the extreme case of an oscillating point potential 
—f35 {x) cos [tot + 77) (see Fig.l). 
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FIG. 1. An tllustration of the system: The incidnet par- 
tides ' energy fl ( and the oscillating frequency uj ) determines 
whether the incident partides will be elevated to htgher energy 
States (dark circles - most cases) or will remain in their initial 
state (open circles - specific cases). 

In terms of the Schròdinger equation, the dynamics 
can be expressed by 

- ip" - (35 {x) cos [ujt + ri) i) + V [x) i/j = (1) 

where we adopt the units 2m ~ l,h = 1, and the no- 
tations ip" = d'allí /dx'^, tp = dip/dt , and V (x) is the 
potential barrier which vanishes quickly for |a;| —> 00. 

The time-dependent solution can be written as a dis- 
crete Fourier transform 

^{x<0)^ ^+e-"* + J2 (2) 

n 



where (p^ are the solutions of the stationary state 
Schròdinger equation, which does not include the oscil- 
lating term, i.e. 

- ^t" + [V{x) - ^] = (4) 

The homogeneous solutions describe waves that 
propagate from ~oo to +00, while íp~ describes the 
waves that are incoming from +00 and outgoing to 
-00. Thus, ip+ ~* ^^giV^x-iut (fQj. X ^ 00) while 

r^e-^^''-''^^ (for x -00) , and |t„|^ is the 
probability of penetrating the barrier. 

By taking care of the matching conditions of the solu- 
tions in eqs. H and ^ at a: = 0, we easily obtain 

2s„Xri + P {Sn-1 + Sn+l) = 2xoV'í'5nO· (5) 

When using the foUowing notations 

(Pn (fin 

and 

s„ = e^^V+í», (7) 

where ip^ = (^o+no. (^^ = 0)- ft' = 5'^í1+„t^/52;U=o, 
(5,iO is the Kronecker delta and í„ = ísi+ni^- 

This difference equation can readily be solved. Notice 
that 

Xn = l/5« (0) (8) 
where gn ix) is the Green function of the equation 

Thus, in the case of a perfectly symmetric rectangular 
barrier, comes directly from |9| 

g-a (0) = - coth [p„L + i arctan (fc„/p„)] /2p„ (9) 

where = v^T+na;, p„ = \JV — k^, 2L is the bar- 
rier length and finally V is its potential height. 

In Fig.2 we present the exact numerical solution of 
eq.|] in the case of a rectangular barrier (in the figure the 
absolute value of s„ is presented), which is the spectral 
solution of eqj^ at a; = 0. This figures illustrates the 
solution's sensitivity to the incoming partides' energy: a 
2% change in Q causes a severe reduction in activation. 

To formulate an analytical expression for this solution, 
we take advantage of the fact that the perturbations are 
strong, i.e., we can assume that ^ F — ^ w. More- 
over, although in the numerical analysis we used the ex- 
act form of the Green function (cq.^) , in the case of an 
opaque barrier, the approximation gn (0) ~ — (2p„) ^ 
may be used with great accuracy. 

In this strong perturbation amplitude and low- 
frequency regime, the difference equation (eq.^ can be 
transformed to a differential equation. By using the def- 
initions 
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G (n) = l3s„/xo 



(10) 



and c (71) = l + Xn//3, one easily obtains the difFcrcntial 
equation 



— Gin) + c{n)G{n)^Sin) 



(11) 



where 5 (n) is the Dirac delta function. 
Hence, one can regard G (n) as having a Green function 
properties (including the singularity at n = 0). 




FIG. 2. A plot of the transmission coefficient |s„| (defined 
in Eq^ as a function of the transmitted partides frequency 
(the activated energy) Q. + nu. The sòlid line represents the 
case Q/V = 0.625 and the dotted one represents the case 
n/V = 0.6125. The other parameters are u/V = 0.0075, 
VVL = 10.75 and /3L = 9.35. 

When the Green function G (n) is known the trans- 
mitted sohition of cqj^ foUows directly froni cqs.||,^ and 



(x > 0) 



f3 



(12) 



In particular, the scattered wave function at a; = is 
proportional to the Fourier transform of the Green func- 
tion. 

Since V — Cl LO, one can approximate 



c~ í~2p/l3 + Lün/ [Pp), 



(13) 



where p = y/V — ü. 

Then, we can define for convenience the variable 



1/3 



and the Green function is then 



(14) 



" Ai(-0 [Ai(-eo)+*Bi(-eo)] 

Ai(-eo)[Ai(-o + íBi(-e)] 



for 
for 



(15) 



where ^0 = ^ (" = 0) = (1 - 2p/f3) {^p/ujy^f^^ and 
Ai and Bi are the Airy functions of the first and second 
kind, respectively. 

For freqüències which are lower than f2 (negative n' s) 
the amplitude oscillates like a simple Airy function (see 
Fig.2 for a typical plot of |s„| , which is related to the 
Green function by eq.p^) 

\Gm^^7r{2(3p/u;f'^,'/'A\'{-0 for n<0 

(16) 

since |Ai(-^o) +íBi(-^o)|^ - n-^^Q^^^- That ex- 
plains the insensitivity of the amplitude of G (for a spe- 
cific n < or ^ < ^0 ) to small variations in the incoming 
energy Q (see, for example, Fig.2). However, for a spe- 
cific incoming energy fi, the amplitude of G oscillates 
with respect to the transmitted eneTgies (i.e., íl + ncj). In 
this regime (i.e., > <^o) 

|G(e)|'^^(2/3pM'/^r'/'Ai2(-eo) for n>0 

(17) 



which means that for an incident partides' energy íí, 
the amplitude of the Green function has a very mild de- 
pendence on the transmitted partides' energies (i.e., on 
n), while it is strongly dependent on the incident parti- 
des' energies íl. This can explain the result presented in 
Fig.2, where a two percent change in the incoming par- 
tides' energy was enough to frustrate the activation to 
higher energies. 

It is clear from eq.|ï^ and from the periodicity of the 
Airy function, that the probability of an incident parti- 
cle being activated to higher energy is very sensitive to 
its initial energy. This sensitivity is manifested in the 
foUowing calculation of the mean activation energy: 



Xofo 
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^n+nu. (0) 



(18) 



A plot of the mean activation energy (ílact) as a func- 
tion of the incident one (f2) is shown in Fig.3 (an exact 
numèric solution). 

The opaquencss of the barrier is the cause for the 
sharp changes in Qact- Since the tunneling coefficient 
(0) 



nu 



{x > L) /(^^ 



+ 



exp i— y/V — íl — nuL) is 
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exponentially small for non-activated partides, there is a 
great advantage for a particle to be activated to higher 
energies + neu ~ V. However, since G (n) is an os- 
cillating function of the incoming energy O, there are 
some energies Vlm for which G {{V — Í7„i) /lo) vanishes. 
In these cases, according to eq. |ï^, not only is the tran- 
sition to energy V forbidden, but the transition to all the 
other higher energies (which corrcspond to n > 0) is, as 
well. 

Therefore, the partides must tunnel out with energy 
which is very close to their initial one (i.e., ííact ~ 
That explains the source of these oscillations and the 
sharp transitions between the activated energies ~ V 
and the non-activated ones ~ íi. The transitions occur 
when the Airy function of eq.|l7| vanishes or, more ac- 
curately, when Ai^ (— Co) — exp[—y/V — i^mL). Taking 
the low-frequency limit (tj — > 0), the Airy function can 
be expressed by simple trigonometric functions to obtain 

V J 4 - 



the transition criterion cos^ 



exp {-pmL) where pm = " 
coming energies for which ílact 
tion occurs, are approximately 



Çljyi- Therefore, the in- 
^ n, and thus no activa- 



V - /3V4 + 



í 3 
-jJUJ I TO + - I TT 



2/3 



(19) 



and the spectral width of these regions is exponentially 
small r ~ exp(— /3L/2). In Fig.3 the non-activated en- 
ergies íim, for which ÍIqcí — are the minima of the 
plot while the màxima correspond to full activation (i.e., 

ílact ^ V). 
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FIG. 3. Two characteristics plots (for arbitrary parame- 
ters) of the mean activation energy (ílact) as a function of 
the incident partides' energy Q, (see Eq\Ld^. The two curves 
represents the same system except for uj, which is five times 
smaller in the dotted curve 

It should be emphasized here that are always much 
larger than the oscillations' frequency lo, and therefore 
the resonances pattern are not related to the oscillator 



quanta (see, for example, ref. pO[|). 

This effect occurs only when the perturbations are 
strong; otherwise (i.e., (3^ <^ y — íl), the solution of 
eq. H is reduced to Sn ^ Iln'^i P/'^Xn', which is always 
(for large n) an exponentially small quantity. In Fig.3 
one can see that when the requirement (3^ ^ V — íl does 
not hold the effect vanishes. 

Moreover, when the oscillations rate decreases (i.e., oj 
decreases) the spectral difference between two successive 
valleys shrinkage, i.e., Aflm = ^m — ^m-i also decreases. 
It is therefore clear, that when lo is small enough so that 
Aíím < r the effect vanishes and elevation to higher 
energies is depressed. In Fig.3 this is shown in the dotted 
line. 

Hence, this effect cannot be anticipated neither in the 
adiabatic regime nor in the weak perturbation approach. 
In both regimes the effect vanishes. 

In general, the exact geometry and shape of the bar- 
rier and the impurity are of no essential importance (see, 
e.g., generalization in ref. |^). One can always (for any 
given geometry) control the activation energy and the 
transmission by adapting the external field frequency. 

To summarize, tunneling transmission through an 
opaque barrier with an oscillating section was investi- 
gated. It was shown that in the strong perturbation (and 
non adiabatic) regime a new sdection-rule appears. Not 
all the incident partides are activated, as could be an- 
ticipated in a tunneling process. Although in most cases 
the incident partides are elevated to higher energy states, 
for some incident partides' energy activation is depressed 
and the partides remain approximately in their initial 
states. The spectral width of these energy domains is 
exponentially small. 

I would like to thank Prof. Mark Azbel for cnlightening 
discussions. 
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